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Abstract 



Two-jet cross sections in deep inelastic scattering at HERA are calculated 
in next-to-leading order. The importance of higher order corrections and 
recombination scheme dependencies is studied for various jet algorithms. Some 
implications for the determination of as(^^), the determination of the gluon 
density and the associated forward jet production in the low x regime at 
HERA are briefly discussed. 
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I. INTRODUCTION 



Deep inelastic scattering (DIS) at HERA is a copious source of multi-jet events. Typical 
two-jet cross sections^ are in the 100 pb to few nb range and thus provide sufficiently high 
statistics for precision QCD tests [0]. Clearly, next-to-leading order (NLO) QCD corrections 
are mandatory on the theoretical side for such tests. Full NLO corrections for one and 
two-jet production cross sections and distributions are now available and implemented in 
the ep — >■ n jets event generator MEPJET, which allows to analyze arbitrary jet definition 
schemes and general cuts in terms of parton 4- momenta 0. A variety of topics can be 
studied with these tools. They include 

• The determination of over a wide range of scales: The dijet cross section is 
proportional to as^i^n) at leading order (LO), thus suggesting a direct measurement of 
the strong coupling constant. However, the LO calculation leaves the renormalization 
scale hr undetermined. The NLO corrections substantially reduce the renormalization 
and factorization scale dependencies which are present in the LO calculations and 
thus reliable cross section predictions in terms of as{mz) (for a given set of parton 
distributions) are made possible. 

• The measurement of the gluon density in the proton (via 7 (7 — >■ qq) : The gluon density 
can only be indirectly constrained by an analysis of the structure function F2 at HERA 
[0]. The boson gluon fusion subprocess dominates the two jet cross section at low x 
and allows for a more direct measurement of the gluon density in this regime. A 
first LO experimental analysis has been presented in NLO corrections reduce the 
factorization scale dependence in the LO calculation (due to the intital state collinear 
factorization, which introduces a mixture of the quark and gluon densities according 
to the Altarelli-Parisi evolution) and thus reliable cross section predictions in terms of 
the scale dependent parton distributions are made possible. 

• The study of internal jet structure: NLO corrections in jet physics imply that a jet (in 
a given jet definition scheme) may consist of two partons. Thus first sensitivity to the 
internal jet structure is obtained, like dependence on the cone size or on recombination 
prescriptions. These studies are also important for reliable QCD studies (such as the 
Us or gluon density determinations). The recombination dependence is only simulated 
at tree level in the NLO calculation and thus the dependence of the cross section on 
the recombination scheme is subject to potentially large higher oder corrections. 

• Associated forward jet production in the low x regime as a signal of BFKL dynamics: 
BFKL evolution |]^ leads to a larger cross section for events with a measured forward 
jet (in the proton direction) with transverse momentum p^''(j) close to Q than the 
DGLAP evolution. Clearly, next-to-leading order QCD corrections for fixed order 
QCD, with Altarelli-Parisi (DGLAP) evolution, are mandatory on the theoretical side 
in order to establish a signal for BFKL evolution in the data. 



^In the following the jet due to the beam remnant is not included in the number of jets. 
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• The determination of the polarized gluon structure function (via 7(7 qq) in polar- 
ized electron on polarized proton scattering 0: The measurement of the polarized 
parton densities and in particular the polarized gluon density from dijet production 
at a polarized electron proton machine at HERA energies would allow to discriminate 
between the different pictures of the proton spin underlying these parametrizations. 
We will analyze these effects in a subsequent publication. 

Some theoretical aspects related to these studies are discussed in the following. 



II. TECHNICAL MATTERS AND JET ALGORITHMS 

The goal of a versatile NLO calculation is to allow for an easy implementation of an arbi- 
trary jet algorithm or to impose any kinematical resolution and acceptance cuts on the final 
state particles. This is best achieved by performing all hard phase space integrals numeri- 
cally, with a Monte Carlo integration technique. This approach also allows an investigation 
of the recombination scheme dependence of the NLO jet cross sections. For dijet production 
at HERA such a NLO Monte Carlo program is MEPJET. 

The basic features of the calculation are described in p[ and we repeat only some of 
them here. In Born approximation, the subprocesses 7* + — >■ Q + g, 7* + ^ — ^ (1 + 9, 
and "y* + g q + q contribute to the two-jet cross section. At 0(a;^) the real emission 
corrections involve 7* + ? 1 + 9 + 9, 1* + 1 1 + Q + 1* + 9 — * 1 + Q + 9 and 
analogous anti-quark initiated processes. The corresponding cross sections are calculated by 
numerically evaluating the tree level helicity amplitudes as given in Ref. p[. The tree level 
matrix elements are numerically checked against the matrix elements in They need 

to be integrated over the entire phase space, including the unresolved regions, where only 
two jets are reconstructed according to a given jet definition scheme. In order to isolate 
the infrared as well as coUinear divergencies associated with these unresolved regions the 
resolution parameter s^m is introduced. This s^m technique has already been successfully 
applied to next-to-leading order calculations of jet cross sections in e~^e~ annihilation and in 



hadronic collisions |Tl|Jl^ . Soft and collinear approximations are used in the region where at 
least one pair of partons, including initial ones, has Sij = 2pi -pj < Smin and the soft and/or 
collinear final state parton is integrated over analytically. Factorizing the collinear initial 
state divergencies into the bare parton distribution functions and adding this soft+collinear 
part to the O(a^) virtual contributions for the 7* + g —>■ q + 9, 1* + q — > 9 + 9, and 
'J* + 9 ^ q + q subprocesses gives a finite result for, effectively, 2-parton final states. In 
general this 2-parton contribution is negative and grows logarithmically in magnitude as 
Smin is decreased. This logarithmic growth is exactly cancelled by the increase in the 3 
parton cross section, once Smin is small enough for the approximations to be valid. 

As mentioned before the collinear initial state divergencies are factorized into the bare 
parton densities introducing a dependence on the factorization scale fip- In order to handle 
these singularities we follow Ref. 0] and use the technique of universal "crossing functions" . 

The integration over the 3-parton phase space with Sij > Smin is done by Monte-Carlo 
techniques (without using any approximations). In general, Smin has to be chosen fairly small 
(below < 0.1 GeV^). Therefore, the effective 2-parton final state (soft+collinear+virtual 
part) does not depend on the recombination scheme. The essential benefit of the Monte 
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Carlo approach in MEPJET is that all hard phase space integrals over the region Sij > 
Smin are performed numerically. Since, at each phase space point, the parton 4-momenta 
are available, the program is flexible enough to implement arbitrary jet algorithms and 
kinematical resolution and acceptance cuts. 

For the numerical studies below, the standard set of parton distribution functions is 



MRS set D-' [|T3[. We employ the two loop formula for the strong coupling constant with 



~ '^^^ MeV, which is the value from the parton distribution functions. The value of is 
matched at the thresholds = nig and the number of flavors is flxed to = 5 throughout, 
i.e. gluons are allowed to split into flve flavors of massless quarks. Unless stated otherwise, 
the renormalization scale and the factorization scale are set to /x/j = = 1/2 J2i Pt(0) 
where p^{i) denotes the magnitude of the transverse momentum of parton i in the Breit 
frame. A running QED fine structure constant aiQ"^) is used. The lepton and hadron 
beam energies are 27.5 and 820 GeV, respectively. A minimal set of kinematical cuts is 
imposed on the initial virtual photon and on the final state electron and jets. We require 
40 GeV^ < < 2500 GeV^, 0.04 < y < 1, an energy cut of E{e') > 10 GeV on the scattered 
electron, and a cut on the pseudo-rapidity r] = — lntan(^^/2) of the scattered lepton and jets 
of |?7| < 3.5. In addition jets must have transverse momenta of at least 2 GeV in both the 
lab and the Breit frame. 

Within these general cuts four different jet definition schemes are considered for which 
we have chosen parameters such as to give similar LO cross sections. Note however, that 
the phase space region for the accepted dijet events depends on the jetalgorithm and thus 
somewhat different event sets would be considered in an actual experiment. Unless stated 
otherwise, and for all jet algorithms, we use the i?-scheme to recombine partons, i.e. the 
cluster momentum is taken as pi +Pj, the sum of the 4-momenta of partons i and j, if these 
are unresolved according to a given jet definition scheme. 

1) W-scheme: 

In the ly-scheme the invariant mass squared, Sij = {pi is calculated for each 

pair of final state particles (including the proton remnant) [Q. If the pair with the 
smallest invariant mass squared is below ycutW^, the pair is clustered according to a 
recombination scheme. This process is repeated until all invariant masses are above 
VcutW^. The resolution parameter i/cut is fixed to 0.02. 

2) JADE-scheme: 

The experimental analyses in [|I| are based on a variant of the H^-scheme, the "JADE" 
algorithm . It is obtained from the py-scheme by replacing the invariant definition 



^ij = iPi + Pi)^ by Mfj = 2EiEj{l — cos 9ij), where all quantities are defined in the 
laboratory frame. Neglecting the explicit mass terms p"^ and p'j in the definition of 
Mfj causes substantial differences in jet cross sections between the W and the JADE 
scheme. 

3) cone scheme: 

In the cone algorithm (which is defined in the laboratory frame) the distance AR = 



[At])^ + (A0)2 between two partons decides whether they should be recombined to a 
single jet. Here the variables are the pseudo-rapidity rj and the azimuthal angle 0. We 
recombine partons with AR < 1. Furthermore, a cut on the jet transverse momenta 
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of PtU) > 5 GeV in the lab frame is imposed in addition to the 2 GeV Breit frame 
cut. 



4) kT scheme: 

For the fc^ algorithm (which is implemented in the Breit frame), we follow the descrip- 
tion introduced in Ref. |jl5|. The hard scattering scale is fixed to 40 GeV^ and 
Vcut = 1 is the resolution parameter for resolving the macro-jets. In addition, jets are 
required to have a minimal transverse momentum of 5 GeV in the Breit frame. 
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FIG. 1. Dependence of the inclusive two-jet cross section in the fey, cone, JADE, and the 
IF-scheme on SjYiiyi^ the two-parton resolution pctrameter. Pctrtons are recombined in the -E'-scheme. 
Error bars represent statistical errors of the Monte Carlo program. For the fairly soft jet definition 
criteria described in the text, Smin independence is achieved for Smin < 0.1 GeV^ 

A powerful test of the numerical program is the s^m independence of the NLO two jet 
cross sections for all jet algorithms. Fig. 1 shows the inclusive dijet cross section as a function 
of Smin for the four jet algorithms. As mentioned before, s^m is an arbitrary theoretical 
parameter and any measurable quantity should not depend on it. One observes that for 
values smaller than 0.1 GeV^ the results are indeed independent of Smin- The strong Smin 
dependence of the NLO cross sections for larger values shows that the soft and collinear 
approximations used in the phase space region Sij < Smin are no longer valid, i.e. terms of 
0{smin) and 0{smin In Smin) bccome important. In general, one wants to choose Smin as large 
as possible to avoid large cancellations between the virtual-|-collinear-|-soft part [sij < Smin) 
and the hard part of the phase space {sij > Smin)- Note that factor 10 cancellations occur 
between the effective 2-parton and 3-parton final states at the lowest Smin values in Fig. 1 
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and hence very high Monte Carlo statistics is required for these points. Smin independence 
is achieved at and below Sm,„ = 0.1 GeV^ and we choose this value for our further studies. 



III. DUET CROSS SECTIONS IN NLO 

A. K-Factors and Recombination Scheme Dependence 

Table |I| shows the importance of higher order corrections and recombination scheme 
dependencies [|16| of the two jet cross sections for the four jet algorithms. While the higher 
order corrections and recombination scheme dependencies in the cone and /ct schemes are 
small, very large corrections appear in the I^-scheme. In addition, the large effective K- 
factor (defined as = <Jnlo/o'lo) of 2.04 (2.02) for the two-jet inclusive (exclusive) cross 
section in the ly-scheme depends strongly on the recombination scheme which is used in the 
clustering algorithm. Such large dependencies are subject to potentially large higher order 
uncertainties, since the recombination dependence is only simulated at tree level in the NLO 
calculation. 

The large corrections and recombination scheme dependencies in particular in the W 
scheme can partly be traced to large single jet masses (compared to their energy in the 
parton center of mass frame). As has been shown in 0, 50 % of the events in the NLO 
cross section for the W scheme (with the E recombination scheme) have a massive jet with 
m/E > 0.44, while substantially smaller values are found in the other jet schemes. The very 
large median value of m/E in the VT-scheme implies that at NLO we are dealing with very 
different types of jets than at LO, and this difference accounts for the large i^-factor. 

In the JADE-algorithm the i^-factor is reduced from 1.48 in the i?-scheme to 1.36 and 
1.24 in the EO and P-scheme^. For the cone (/ct) scheme this recombination scheme de- 
pendence is reduced to the 3% (10%) level. 



TABLE I. Two-jet cross sections in DIS at HERA. Results are given at LO and NLO for the 
four jet definition schemes and acceptance cuts described in the text. The 2-jet inclusive cross 
section at NLO is given for three different recombination schemes. 





2-jet 


2-jet exchisive 


2-jet inclusive 


2-jet inclusive 


2-jet inclusive 




LO 


NLO (E) 


NLO (E) 


NLO (£^0) 


NLO (P) 


cone 


1107 pb 


1047 pb 


1203 pb 


1232 pb 


1208 pb 


kT 


1067 pb 


946 pb 


1038 pb 


1014 pb 


944 pb 


W 


1020 pb 


2061 pb 


2082 pb 


1438 pb 


1315 pb 


JADE 


1020 pb 


1473 pb 


1507 pb 


1387 pb 


1265 pb 



^The NLO two jet cross sections for the W or the JADE scheme in Table I disagree with previous 
calculations 



The DIS JET program [10|, for example, gives a i^T-factor very close to unity for 
a phase space region which is very similar. 
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FIG. 2. Dependence of a) the two-jet inclusive and b) the two-jet exclusive cross section in 
the cone scheme on the renormalization and factorization scale factor ^. The solid curves are for 
/^?? = /^F = ? Pt(0)^' while for the dashed curves only £,r = C varied but = 1/4 is fixed. 
Choosing the photon virtuality as the basic scale yields the dotted curves, which correspond to 
/^fl — l^^F ~ ^ Q"^- Results are shown for the LO (lower curves) and NLO calculations. 



B. Scale Dependence 



As mentioned before, the NLO corrections substantially reduce the the renormalization 
and factorization scale dependence which is present in the LO calculations and thus reliable 
cross section predictions in terms of asijriz) are made possible. The scale dependence for 
dijet cross sections in the cone scheme is shown in Fig. 2. We have considered scales related 
to the scalar sum of the parton transverse momenta in the Breit frame, Yl,i Pt{'^)i ^"^^ the 
virtuality of the incident photon. In Fig. |^ the dependence of the two-jet cross section, in 
the cone scheme, on the renormalization and factorization scale factors and E,f is shown. 
For scales related to Pr(^) ^^ey are defined via 

i i 

For the two-jet inclusive cross section of Fig. the LO variation by a factor L43 is reduced 
to a 10% variation at NLO when both scales are varied simultaneously over the plotted 
range (solid curves). However, neither the LO nor the NLO curves show an extremum. The 
uncertainty from the variation of both scales for the NLO two-jet exclusive cross section in 
Fig. 2b (solid curves) is reduced to 5%. Furthermore, the two-jet exclusive cross section now 
has a maximum and is equal to the LO cross section for ^ = 0.5. Also shown is the ^ = 
dependence of LO and NLO cross sections at fixed = 1/4 (dashed curves). In this case 
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a maximum appears in the NLO inclusive and exclusive cross sections. However, the scale 
variation is stronger than in the ^ = = C,f case. 

An alternative scale choice might he ^\ = = ^Q^ . The resulting ^ dependence is 
shown as the dotted lines for both the LO and NLO calculations. At LO the two scale choices 
give quahtatively similar results. However, with f^R = fij? = C, , the scale dependence does 
not markedly improve at NLO. We therefore use the jet transverse momenta in the Breit 
frame to set the scale and fix = (^j? = 1/4 in Eq. (|I|) for the following numerical results. 
A careful study of the scale dependence and the choice of the scale in the dijet cross section 
is needed in order to extract a reliable value for as(M^). 



IV. GLUON DENSITY DETERMINATION 

HERA opens a new window to measure the proton structure functions, in particular 
the gluon distribution, in a completely new kinematic region. The accessible range in the 
Bj or ken-scaling variable x can be extended considerably towards low x compared to previous 
fixed target experiments. Dijet production in DIS at HERA in principle allows for a more 
direct measurement of the gluon density in the proton (via 7(7 qq) than an analysis of 
the structure function F2. 

For these studies we use the cone scheme as defined in section II. The range is 
lowered to 5 < < 2500 GeV^ and the cut on the jet transverse momenta in the Breit 
frame is increased to 5 GeV. The LO (NLO) results are based on the LO (NLO) parton 



distributions from GRV [0 together with the one-loop (two-loop) formula for the strong 
coupling constant. With these parameters, one obtains 2890 pb (2846 pb) for the LO (NLO) 
two jet exclusive cross section. 

In order to investigate the feasibility of the parton density determination at low x. Fig. 3a 
shows the Bjorken x distribution of the two jet exclusive cross section in the cone scheme. 
The gluon initiated subprocess clearly dominates the Compton process for small x in the LO 
predictions. The effective ii"-factor close to unity for the total exclusive dijet cross section 
is a consequence of compensating effects in the low x {K > 1) and high x {K <1) regime. 

For the isolation of parton structure functions we are interested in the fractional mo- 
mentum Xi of incoming parton i {i = q,g), however, and in dijet production x and Xi differ 
substantially. Denoting as sjj the invariant mass squared of the produced dijet system, and 
considering two-jet exclusive events only, the two are related by 

x, = x(l+'-^] (2) 



The Sjj distribution of Fig. 3b exhibits rather large NLO corrections as well. The invariant 
mass squared of the two jets is considerable larger at NLO than at LO (the mean value of 
Sjj rising to 570 GeV^ at NLO from 470 GeV^ at LO). 

The NLO corrections to the x and Sjj distributions have a compensating effect on the 
Xi distribution in Fig. 3c, which shows very similar shapes at LO and NLO. At LO a direct 
determination of the gluon density is possible from this distribution, after substraction of 
the calculated Compton subprocess. This simple picture is modified in NLO, however, and 
the effects of Altarelli-Parisi splitting and low pt partons need to be taken into account more 
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FIG. 3. a) Dependence of the exclusive two-jet cross section in the cone scheme on Bjorken 
X for the quark and gluon initiated subprocesses and for the sum. Both LO (dashed) and NLO 
(soUd) results are shown; b) Dijet invariant mass distribution in LO (dashed) and in NLO (solid); 
c) Same as a) for the Xj distribution, Xi representing the momentum fraction of the incident parton 
at LO. 
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TABLE 11. Forward jet cross sections in [pb] in DIS at HERA. 





W i tli 


\K7^ "I" h on "I" 

W 1 UilU U. h 


1 C-ldLl V C 




forward jet 


forward jet 


phase space 


1 jet (LO) 


pb 


9026 pb 


0% 


2 jet (LO) 


19.3 pb 


2219 pb 


0.87% 


2 jet (NLO) 


68 pb 


2604 pb 


2.61% 


3 jet (LO) 


30.1 pb 


450 pb 


6.7% 



carefully to determine the structure functions at a well defined factorization scale fip. We 
will further investigate this problem in a subsequent publication. One method to determine 
the gluon density in NLO is presented in |jT9[ . 



V. FORWARD JET PRODUCTION IN THE LOW X REGIME 

Deep inelastic scattering with a measured forward jet with relatively large momentum 
fraction Xj^t (in the proton direction) and PT'^^{j) ~ is expected to provide sensitive 
information about the BFKL dynamics at low x [^,^. In this region there is not much 
phase space for DGLAP evolution with transverse momentum ordering, whereas large ef- 
fects are expected for BFKL evolution in x. In particular, BFKL evolution is expected to 



substantially enhance cross sections in the region x « Xjet p0| , |21|| . In order to extract 
information on the ln(l/x) BFKL evolution, one needs to show that cross section results 
based on fixed order QCD with DGLAP evolution are not sufficient to describe the data. 
Clearly, next-to-leading order QCD corrections to the DGLAP predictions are needed to 
make this comparison between experiment and theory. 

In Table we show numerical results for the multi jet cross sections with (or without) a 
forward jet. The LO (NLO) results are based on the LO (NLO) parton distributions from 



GRV |jT8|] together with the one-loop (two-loop) formula for the strong coupling constant. 
Kinematical cuts are imposed to closely model the HI event selection |^ . More specifically, 
we require Q"^ > 8 GeV^ , x < 0.004, 0.1 < y < 1, an energy cut of E{e') > 11 GeV on 
the scattered electron, and a cut on the pseudo-rapidity rj = — lntan(^/2) of the scattered 
lepton of —2.868 < ■r]{e') < —1.735 (corresponding to 160° < 9{l') < 173.5°). Jets are defined 
in the cone scheme (in the laboratory frame) with AR = 1 and |?7(j)| < 3.5. We require a 
forward jet with Xjet = Pzij)/Ep > 0.05, E{j) > 25 GeV, 0.5 < PTij)/Q^ < 4, and a cut on 
the pseudo-rapidity of 1.735 < ri{j) < 2.9 (corresponding to 6.3° < 9{j) < 20°). In addition 
all jets must have transverse momenta of at least 4 GeV in the lab frame and 2 GeV in the 
Breit frame. 

The cross sections of Table || demonstrate first of all that the requirement of a forward 
jet with large longitudinal momentum fraction {xjet > 0.05) and restricted transverse mo- 
mentum (0.5 < PT{j)/Q'^ < 4) severely restricts the available phase space, in particular for 
low jet multiplicities. The 1-jet exclusive cross section vanishes at LO, due to the contradict- 
ing X < 0.004 and Xjet > 0.05 requirements. For x << Xjet, a high invariant mass hadronic 
system must be produced by the photon-parton collision and this condition translates into 



10 



2E{j)mT ^ S^,parton ~ " l) » ' (3) 

where and y are the transverse mass and rapidity of the partonic recoil system, respec- 
tively. Thus a recoil system with substantial transverse momentum and/or invariant mass 
must be produced and this condition favors recoil systems composed out of at least two 
additional energetic partons. 

As a result one finds very large fixed order perturbative QCD corrections (compare 2 
jet LO and NLO results with a forward jet in Table In addition, the LO {0{a^)) 3-jet 
cross section is larger than the LO {0{as)) 2-jet cross section. Thus, the forward jet cross 
sections in Table H are dominated by the {0{al)) matrix elements. The effects of BFKL 
evolution must be seen and isolated on top of these fixed order QCD effects. We will analyze 
these effects in a subsequent publication. 

VI. CONCLUSIONS 

The calculation of NLO perturbative QCD corrections has received an enormous boost 
with the advent of full NLO Monte Carlo programs [Tl|JT^j2^ . For dijet production at HERA 
the NLO Monte Carlo program MEPJET allows to study jet cross sections for arbitrary 
jetalgorithms. Internal jet structure, parton/hadron recombination effects, and the effects 
of arbitrary acceptance cuts can now be simulated at the full 0{a'^) level. We found large 
NLO effects for some jet definition schemes (in particular the ly-scheme) and cone and kx 
schemes appear better suited for precision QCD tests. 

The extraction of gluon distribution functions is now supported by a fully versatile NLO 
program. Preliminary studies show that large NLO corrections are present in the Bjorken x 
distribution for dijet events, while these effects are mitigated in the reconstructed Feynman 
X (xi) distribution, thus aiding the reliable extraction of g{xi,nj^). 

For the study of BFKL evolution by considering events with a forward "Mueller"-jet 
very large QCD corrections are found at (9(a^). These fixed order effects form an important 
background to the observation of BFKL evolution at HERA. They can now be studied 
systematically and for arbitrary jetalgorithms. 
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